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О  хМОМЕНТЕ  КОЛИЧЕСТВА  ДВИЖЕНИЯ  СОЛНЦА 

Сделан  расчет  распределения  плотности,  массы  н  момента  количе- 
ства движения  для  отдельных  слоев  Солнца  и  части  его,  заключенной 
внутри  заданной  доли  радиуса;  при  этом  была  принята  модель  Шварц- 
шильда,  дающая  малые  отличия  в  отношении  величины  радиуса  инерции 
Солнца  и  его  момента  количества  движения  Ь  по  сравнению  с  моделями, 
предложенными  после.  Согласно  этой  модели,  99%  массы  Солнца  заклю- 
чено в  сфере,  с  радиусом  г=0,60  Е^^  и  90%  —  в  сфере  с  г=0,41  К^. 
90%  момента  количества  движения  Солнца  (считая  его  вращающимся 
как  твердое  тело)  заключено  в  сфере  радиусом  г=0,55  Е^.  Если  рассма- 
тривать сферические  слои  внутри  Солнца  одинаковой  толщины,  то  наи- 
оольшей  массой  будет  обладать  слой  с  радиусом  около  0,2  ^?0,  а  наи- 
большим  моментом  количества  движения  —  слой  с  г=0,3  Е^. 

Сравнивая  момент  количества  движения  Солнца,  вращающегося  как 
твердое  тело,  с  современным  моментом  количества  движения  всех  планет, 
получаем,  что  момент  количества  движения  Солнца  в  этом  предположе- 
нии составляет  только  0,45%  общего  момента  количества  движения 
планет.  Если  значительная  часть  протопланетного  материала  падала 
на  Солнце,  то  она  должна  была  сообщить  Солнцу  больший  момент  коли- 
чества движения,  нежели  наблюдается  сейчас.  Это  обстоятельство  яв- 
ляется одним  из  решающих  затруднений  теории  Лапласа  и  должно  учи- 
тываться космогоническими  гипотезами. 

В  статье  рассматривается  возможность  наличия  большего  момента 
г;оличества  движения  у  Солнца  при  вращении  его  внутренних  частей 
1-  большей  угловой  скоростью,  нежели  поверхностные  слои.  Расчет  носит 
характер  приближенной  оценки  максимального  возможного  момента 
количества  движения.  При  этом  использзгются  наблюдательные  данные 
относительно  сжатия  поверхностных  слоев  Солнца,  которое  принимается 
меньшим  10~^,  что  соответствует  разнице  между  экваториальным  и  поляр- 
ным радиусами  Солнца,  меньшей  0','1. 

Рассмотрено  несколько  моделей  строения  Солнца.  Для  модели, 
наиболее  близкой  к  динамически-равновесной,  максимальный  допусти- 
мый момент  количества  движения  Солнца  полз'чается  всего  в  4,3  раза 
Польше,  чем  для  модели  Солнца  Шварцшильда  и  вращающегося  как 
гвердое  тело.  Таким  образом,  даже  в  этом  предельном  случае  момент 
количества  движения  Солнца  составляет  всего  около  2%  общего  момента 
количества  движения   всех   планет  солнечной   системы. 


//.  //.  Парийский 


Этим  исключается  возможность  выпадения  на  Солнце  сколько-нибудь 
значительной  доли  материала,  послужившего  для  образования  планет, 
если  этот  материал  вращался  вокруг  Солнца  (подобно  планетам)  преиму- 
щественно в  одну  сторону  и  если  не  существовало  какого-либо  эффектив- 
ного механизма  потери  Солнцел!  его  вращательного  момента  количества 
движения. 

§  1.  Вращение  Солнца  как  твердого  тела 

Для  ряда  космогонических  гипотез  существенно  знать,  какая 
часть  момента  количества  движения  солнечной  системы  может 
быть  заключена  в  Солнце.  Момент  количества  движения  сол- 
нечной системы,  с  учетом  близости  плоскостей  планетных 
орбит  между  собой  и  малости  моментов  количества  движения 
орбитального  движения  спутников  и  вращательного  движения 
самих  планет  относительно  своих  осей,  может  быть  записан 
в  таком  виде: 

10 

где  ^А.  —  мохмент  количества  движения  планеты  к  (считая 
астероиды  за  одну  планету);  Ь©  —  момент  количества  движе- 
ния Солнца;  Ь  —  момент  количества  движения  всех  планет. 
Если  Солнце  вращается  как  твердое  тело  с  угловой  ско- 
ростью (О,  то  его  момент  количества  движения  будет 

•"о 

^0 = ^0 " = "з  "^'^  1  ^'^^^  ^''^  ^^'  ^^^ 

о 

где  /©  —  момент  инерции  Солнца  относительно  оси  вращения; 
г©  —  радиус  Солнца;  р  (г)  выражает  закон  распределения 
плотностей  внутри  Солнца. 

Обозначив  массу  Солнца  через  М©  и  радиус  инерции  Солнца 
относительно  оси  вращения  через  /©,  мы  можем  написать 

■^©=9 -5-^0''©*'^= ^©4"'  (3) 

где 

'■0 

Г    г*р  (г)  (1г 

-т 


5     1     о  ^    ^  ^ 


3  г^©  ;0  ^  \^©- 

I     г2р(г)г7/- 
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Для  случая  постоянной  плотности,  как  хорошо  известно, 
д=1,00;     4  =  |4;     -^  =  0,632. 

Если  мы  обозначим  через  М  (г),  Ь  (г),  д  (г)  соответственно 
массу,  момент  количества  движения  и  множитель  д,  относя- 
щиеся к  части  Солнца,  заключенной  внутри  радиуса  г,  то 

Ь{г)  =  д{г)^М{г)?^^  (3') 


</г 


^ 

г 

I  г2р  (Г)  аг 


^('•)=1г44 •  (4') 


Как  известно,  различные  модели  внутреннего  строения 
Солнца  дают  сравнительно  мало  различающиеся  распределения 
плотностей  с  центральной  плотностью  около  100  г/см^.  По- 
скольку нас  интересует  лишь  приблизительный  расчет,  мы  не 
будем  вдаваться  в  детали  этого  вопроса,  а  воспользуемся 
последней  моделью  Солнца,  предложенной  Шварцшильдом  [1]. 
Согласно  модели.  Солнце  состоит  из  оболочки,  находящейся 
в  лучистом  равновесии,  и  конвективного  ядра  с  радиусом 
около  0,12  радиуса  Солнца.  Химический  состав  Солнца  для 
этой  модели  был  получен  следуюпщй:  47%  водорода,  41%  ге- 
лия и  12%  ресселовской  смеси  остальных  элементов.  В  этой 
модели  центральная  плотность  равна  111,6  г/см^,  или  79,1 
средней  плотности. 

По  данным,  приведенным  в  работе  Шварцшильда,  и  табли- 
цам структуры  политропного  газового  шара  с  показателем 
политропы  л =3/2,  мы  получили  распределение'  плотностей, 
а  затем  масс  и  моментов  количества  движения  для  отдельных 
слоев.  Результаты  приведены  в  табл.  1  и  на  рис.  1. 

В  первом  столбце  таблицы  приведено  расстояние  от  центра 
Оэлнца  г  в  долях  радиуса  Солнца;  во  втором  —  плотность 
р  (г);  в  третьем  —  масса  Л/ (г),  заключенная  внутри  сферы 
радиуса  г,  в  граммах,  деленная  на  1С;  в  четвертом  —  момент 
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инерции  /  (/•)  массы  М  (г)  относительно  оси  вращения  в  системе 
С08,  деленный  на  10^^;  в  пятом  —  величина  д\  в  шестом  — 
плотность  р^,  принимая  за  единицу  среднюю  плотность  Солнца; 
в  седьмом  и  восьмом  —  величины  М  {г)  и  /  (г),  выраженные 
в  долях  их  значений  для  всего  Солнца. 

При  этой  модели  Шварцшильда  для  Солнца  в  целом: 

д  =  0,1284;      -^=0,2267. 


^^  м  м.цйМр 


Рис.     1.     Распределение  механических  характеристик  внутри  Солнца 
согласно  модели  Шварцишльда 

7—  р(г);       2  —  ДМ(г);  3  —  М,  (г);    <  —  4/ (г);    6  —  1Лг)\ 
7  —  Я  (г);    [/  (г)=9  1  г^  М(г)] 

Для  модели,  рассчитанной  И.  Эпштейном,  учитывающей 
протон-протонную  реакцию  и  реакцию  Бете  с  новыми  значе- 
ниями постоянных,  получается  значительно  отличное  отно- 
сительное содержание  элементов  [2],  но  значение  /©/г©  полу- 
чается почти  таким  же,  как  и  для  модели  Шварцшильда. 
Именно,  по  подсчетам  Л.  Мотца   [3], 

^:=:0,123     и      -^-  =  0,222. 

'"о 

Возможные    изменения    этой   модели   вряд   ли   смогут    сильно 
изменить  значение   радиуса   инерции. 
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Таблица    1 

Распределение   . 

неханических 

характеристик    внутри    Солнца  в  модели 

Шварцшильда 

г 

Р 

ДГ(г)10~" 

I  (г)  •  1 0~*' 

Ч 

Р1 

МЛг) 

Л  (г) 

0,00 

111,6 

0,000 

0,000 

1,00 

79,1 

0,000 

0,0000 

0,05 

105,9 

0,019 

,0,001 

0,99 

75,1 

0,010 

0,0003 

0,10 

90,1 

0,138 

0,026 

0,98 

63,9 

0,070 

0,005 

0,15 

68,0 

0,400 

0,162 

0,94 

48,2 

0,201 

0,083 

0,20 

43,8 

0,754 

0,517 

0,89 

31,0 

0,380 

0,105 

0,25 

25,07 

1,110 

1,097 

0,82 

17,77 

0,559 

0,223 

0,30 

13,33 

1,402 

1,807 

0,74 

9,45 

0,707 

0,367 

0,35 

6,815 

1,616 

2,528 

0,66 

4,830 

0,814 

0,514 

0,40 

3,416 

1,760 

3,176 

0,58 

2,421 

0,887 

0,645 

0,45 

1,695 

1,852 

3,709 

0,51 

1,201 

0,934 

0,753 

0,50 

0,8346 

1,910 

4,121 

0,45 

0,5915 

0,962 

0,837 

0,55 

0,4018 

1,944 

4,422 

0,39 

0,2848 

0,980 

0,898 

0,60 

0,1897 

1,963 

4,628 

0,34 

0,1344 

0,990 

0,940 

0,65 

0,0875 

1,974 

4,763 

0,30 

0,0620 

0,995 

0,967 

0,70 

0,0371 

1,980 

4,845 

0,26 

0,0263 

0,998 

0,984 

0,75 

0,01457 

1,982 

4,890 

0,23 

0,0103 

0,999 

0,993 

0,80 

0,00486 

1,984 

4,911 

0,20 

0,00344 

0,9997 

0,998 

0,85 

0,00129 

1,984 

4,923 

0,18 

0,00091 

1,0000 

1,000 

0,90 

0,00024 

1,984 

4,923 

0,16 

0,00017 

1,0000 

1,000 

0,95 

0,00005 

1,984 

4,923 

0,14 

0,00004 

1,0000 

1,000 

1,00 

0,00000 

1,984 

4,923 

0,13 

0,00000 

1,0000 

1,000 

Считая,  что  Солнце  вращается  как  твердое  тело,  мы  для 
момента    количества     движения     Солнца     получим    величину 

^^  =  7^0)0  =  71/о4о)о  =  1,412.  10*8  С08, 
где  принято 

3/©  =  1,984. 1033  г;     /^  =  0,226  г©;     Г0  =  6,95  •  10^'^  см; 
-о=864етб  =  2,868. 10-0  сек- 

Общий  момент  количества  движения  планет  солнечной 
системы  (если  пренебречь  наклонами  плоскостей  орбит  к  неиз- 
менной плоскости  Лапласа)  равен 


10 
к— 1 


(6) 
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где 


Ь,  =  \/§аМе  т^  \/1  +  т^  \/«л  (^  —  4)  =  ^/ёоЩ  т^  1'^  .       (6') 

Здесь  §  —  постоянная  тяготения;  а  —  большая  полуось  земной 
орбиты;  /Ид.  —  масса  планеты  в  граммах;  /п^.©  —  масса  пла- 
неты в  долях  солнечной;  а\  —  большая  полуось  планетной 
орбиты  в  долях  астрономической  единицы;  е^  —  эксцентриситет 

планетной  орбиты;  Ь\.=\/1  +  /^/^-о  \/<^'а-(1  —  б\)- 

Численные  значения  этих  величин  приведены  в  табл.  2. 

Таблица   2 
Момент  количества  дви^кения  планет  солнечной  системы 


Планета 

а 

е 

т/М©- 10» 

ь' 

^.10-19 
т 

1.-10— 2^ 

Меркурий    .    .      0,387 

0,2056 

0,133 

0,61 

2,71 

0,00007 

Венера  .    . 

I    0,723 

0,0068 

2,45 

0,84 

3,74 

0,00182 

Земля    .    . 

\    1,000 

0,0167 

3,00 

1,00 

4,45 

0,00265 

Марс  .    .    . 

1    1,524 

0,0933 

0,323 

1,23 

5,47 

0,00035 

Астероиды 

1    2,910 

0,1000 

0,299 

1,70 

7,55 

0,00045 

Юпитер    . 

1    5,203 

0,0484 

955 

2,28 

10,15 

1,92236 

Сатурн  .   . 

9,539 

0,0558 

286 

3,08 

13,72 

0,77836 

Уран  .    .    . 

19,191 

0,0471 

43,7 

4,38 

*9,46 

0,16875 

Нептун     . 

30,071 

0,0086 

51,8 

5,48 

24,39 

0,25076 

Плутон     . 

39,518 

0,2468 

0,300 

6,09 

27,10 

0,00161 

Отсюда 


2  =  1343  .  10-^ 


^р  =  312,7.  10^8  С08, 


и  общая  масса  планет 


10 


2  =  3,127  •  10 


т. 


=  2  т;^  =  1,343-  Ю-т, 


Отсюда  получаем 


Л=1 


0,9987; 


50 


Л/©  +  т, 


=  0,00134 


о  моменте  количества  двимсения  Солнца  11 


^ 


^Г^+^-г 


0,00450  «0,5%. 


Момент  количества  движения  Солнца  составляет,  таким 
образом,  менее  0,5%  общего  момента  количества  движения 
солнечной  системы. 


§  2.  Наблюдательные  данные 

Попытаемся  выяснить,  насколько  наблюдаемое  сжатие  сол- 
нечной поверхности  может  дать  сведения  о  вращении  внутрен- 
них частей  Солнца.  Разница  между  экваториальным  и  поляр- 
ным радиусами  Солнца  измерялась  многократно. 

Амброн  и  Шур  по  наблюдениям  с  гелиометром  в  течение 
12  лет  не  обнаружили  разницы  между  экваториальным  и 
полярным   радиусами   с   точностью   до    п:0»007. 

Шевалье   по   наблюдениям   за  5  лет  (1905 — 1910)   получил 

а  — с  =  — 0",19±0",02. 

Хайн  по  фотографическим  наблюдениям,  ослабляя  яркость 
Солнца,  получил 

а  — с  =  — 0",035  ±  0",014. 

Если  учесть  изменение  склонения  Солнца  за  время  экспозиции 
в  3  сек.,  как  это  сделал  Шауб,  то  эта  разница  становится 
меньше  точности  наблюдений. 

По  последней  работе  Мейермана,  заново  обработавшего 
наблюдения  Амброна  и  Шура,  производившиеся  с  гелиометром 
с  1890  по  1902  г.,  средний  радиус  Солнца  для  эпохи  1893,5 
равен  Л0=959",935  + О", 0025  с  колебаниями  порядка  О", 05.  Рас- 
хождения измерений  могут  вызываться  скорее  всего  измене- 
ниями в  строении  хромосферы,  имеющей,  как  известно,  струй-- 
чатый  характер.  Поэтому,  хотя  измерения  разности  эквато- 
риального и  полярного  радиусов  Солнца  производятся  с  чрез- 
вычайно высокой  точностью,  мы  должны  считаться  с  возмож- 
ной разницей  в  строении  полярной  и  экваториальной  областей 
хромосферы,  а  возможно,  и  фотосферы. 

Вальдмайер  считает,  что  разность  между  экваториальным 
и  полярным  радиусами  Солнца  не  может  превосходить  0",01  14]. 
Однако,  учитывая  возможные  изменения  в  строении  хромо- 
сферы, частично  сказывающиеся  и  на  строении  фотосферы,  мы 
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будем  считать,  что  эта  разность  не  превосходит  0",1,  т.  е.  при- 
мем в  10  раз  большую  величину,  нежели  Вальдмайер.  Таким 
образом,  мы  будем  считать,  что  сжатие  Солнца 

^0—        ал       ^10000  •  ^   '' 


§  3.  Влияние  центробежной  силы 

Нас  интересует  вопрос,  не  может  ли  внутри  Солнца  нахо- 
диться сжатое  вращающееся  ядро,  которое  совместно  с  осталь- 
ными слоями  Солнца  давало  бы  внешнюю  уроненную  поверх- 
ность со  сжатием,  меньшим  10~*. 

Прямым  эффектом  потенциала  центробежной  силы  мы  можем 
пренебречь,  так  как  этот  потенциал  должен  относиться  к  той 
угловой  скорости,  которая  имеется  во  внешних,  самых  поверх- 
ностных слоях  Солнца  и  практически  не  может  изменить  его 
уроненную  поверхность. 

Действительно,  пусть  потенциал  тяготения  невращающе- 
гося  Солнца  на  его  поверхности  равен 

где  §  —  постоянная  тяготения,  М  —  масса  Солнца,  г^  —  его 
радиус. 

Если  бы  вращались  только  поверхностные  слои  Солнца, 
содержащие  ничтожно  малую  долю  общей  массы  Солнца,  то 
уроненная  поверхность,  определяющая  форму  Солнца,  имела  бы 
вид: 

V,  {Го  +  Ч)  =  ~^^  +  4-  "'  (^'  +  У")  =  'о . 

где  §Г(,  —  отклонение  уровенной  поверхности  от  сферы.  На 
полюсах  Солнца  эта  уровенная  поверхность  совпадала  бы 
с  уровенной  поверхностью  силы  тяготения  (при  том  же  с^), 
а  на  экваторе  она  проходила  бы  дальше  от  центра  Солнца  на 
величину  ^Гц,   определяемую  из  уравнения 

Но 
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Поэтому  ДЛЯ  8Гц  получим  выражение; 


1 


со^г! 


^'•о  =  ,-^^-^='^|^  =  7,2.10^см  =  7,2км. 


Но  на  Солнце  1"  соответствует  726  км.  Следовательно,  вслед- 
ствие учета  потенциала  центробежной  силы  экваториальный 
радиус  Солнца  изменяется  всего  на  О", 01,  так  что  его  влиянием 
можно  пренебречь. 

Если  бы  внутренние  слои  Солнца  вращались  быстрее,  то, 
независимо  от  закона  их  вращения,  внешняя  уроненная  поверх- 
ность все  равно  определялась  бы  потенциалом  тяготения  вну- 
тренних частей  и  потенциалом  центробежной  силы,  соответ- 
ствующей угловой  скорости  вращения  поверхностных  слоев. 
Поэтому  влияние  этого  потенциала  все  равно  было  бы  мало. 

Внешняя  уроненная  поверхность  Солнца  будет,  таким  обра- 
зом, определяться  практически  только  потенциалом  силы  тяго- 
тения и  зависеть  от  распределения  плотностей  внутри  Солнца 
и  сжатия  поверхностей  одинаковой  плотности.  Влияние  рас- 
пределения скоростей  вращения  внутри  Солнца  на  внешнюкт 
уроненную  поверхность,  таким  образом,  является  не  прямым, 
а  косвенным.  Оно  определяется  сжатием  внутренних  слоев, 
которое  в  свою  очередь  определяется  распределением  угловых 
скоростей. 

§  4.  Однородные  сфероиды,  дающие    эквипотенциальные 

поверхности  силы  тяготения,  практически 

совпадающие  с  поверхностью  Солнца 

Рассмотрим  вначале  вспомогательную  задачу.  Предположим, 
что  внутри  Солнца  находится  сжатый  однородный  эллипсоид 
вращения,  центр  которого  совпадает  с  центром  Солнца  и  эква- 
ториальная плоскость  —  с  экваториальной  плоскостью  Солнца. 
Пусть  большая  полуось  этого  сфероида  равна 

а  =  ка^,  (9) 

а  малая 

6  =  а(1  — с)  =  А:ао(1— с),  (90 

где  А  :^  1,  а^  —  радиус  солнечного  экватора,  с  —  сжатие  вну- 
треннего эллипсоида. 

Определим,  каково  должно  быть  сжатие  этого  эллипсоида  с 
в   зависимости  от  Н,   чтобы  внешняя   уроненная   поверхность 
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СИЛЫ  тяготения  этого  эллипсоида,  касающаяся  поверхности 
Солнца  по  его  экватору,  имела  сжатие  с^,  равное  максималь- 
ному сжатию  поверхности  Солнца  Сц. 

Как  хорошо  известно  [5],  потенциал  тяготения  однородного 
эллипсоида  вращения  с  полуосями  а  и  6  на  внешнюю  точку 
Р  (х,  у,  г)  может  быть  написан  в  виде 

У,  =  К[С{1)-В{А){х'^-\-7/)-Л{\)гЦ,  (10) 

где  X  —  больший  корень   уравнения 

а:2  +  у-     I 


Кроме  того, 


«2  + А       '     62+Х 


=  1.  (10а) 


С  (л)  =  («2  —  6^)  агс  1^1;  (10в) 

В{1)=Цагс1§1-^);  (Юг) 

А{А)  =  1  —  атсЬ§1  (10д) 


Внешняя  потенциальная  поверхность  этого  эллипсоида 
будет  некоторым  сжатым  овалоидом  (строго  говоря,  не  эллип- 
соидом) с  экваториальным  и  полярным  радиусами,  равными 
соответственно  а^  и  Ь^. 

Будем  измерять  «сжатие»  этой  поверхности  величиной 

.'_«х-Ьх      д      ^^^^3^=^^  (11) 


«X 


Рассмотрим  сперва  точку  уровенной  поверхности  на  экваторе 
Солнца.  Для  нее 

х^-\-у^  =  а^     и     2  =  0. 
Тогда  из  (10а) 

и 

л  =  а2(1— А;2).  '  (12) 


Из  (10е) 
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/— -        ^^ П2а) 

У^1-А2е2'  ^        ' 


где  е  —  эксцентриситет  сфероида 

из  (106) 

^:=:2гр^^5^.  ■  (13) 

Из  (10),  (10в)  и  (Юг)  имеем,  полагая  ке  =  1: 

\\  =  Ка1  [^2  агс  1д  /  -  агс  1^  /  +  -:^^^']  ,  (14) 

или,  разлагая  в  ряд  до  пятых  степеней  /  и  учитывая,  что  I 
и  е  одного  порядка: 

Но 

^  =  ц^^^^)-"^  =  ^^^^^^^,,,;  (15) 

/з  =  гз  +  |/^+...;  (15а) 

1'^=:1-^-{- ...  (156) 

Поэтому 

Окончательно  для  точки  на  экваторе: 

для  точки  на  полюсе  Солнца 

Тогда  из  (10а) 

а1{1-с,)  =  а1к'е'  +  1  (17) 
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и  из  (10е) 

И 

г  =А:а2[«2агс1д/-(1-д2(г-агс1д/)],  (18) 


или 


У=Ка1[е[1-{р)-а-с,у[1-1  +  ^-\1^)], 


ИЛИ 


I  \  =  Ка1  \1'  (1  -  с,)-^  -  4  ^М1  -  ^о)"'  -^^Ч^-  СоГ'  + 

откуда 

V,  =  Ка1  [1 1^  (1  -  с,)-^  - 1  I'  (1  -с^)-^]  , 


или  окончательно 


\\  =  ^Ка11Ц\-сГ'[}-  УЧ\-С,)-''].  (19) 

Приравнивая  выражения  (16)  и  (19)  для  потенциалов  на 
полюсе  и  экваторе  Солнца,  получаем  уравнение  для  опреде- 
ления I: 

1  +  4^'  =  (1-^о)-'[1-|^'(1-Со)-^],  (20) 

или 

Отсюда  получаем: 

^  1     ,  1  —  3  ^0  •  ^-'^' 

10"^  5(1  — Со)2 

Сжатие  внутреннего  сфероида  связано  с  его  эксцентриси- 
тетом выражением 

а        2  ^2 
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2  _  10  Со 
3   А2" 

^_  5    Со 
^■^  3  А2  ' 

^=т\/Ь 

^0 

е2 

_10  Со         5 
3   А2  ■—  3 

е1 

Л2 

НО 


(22) 

(22а) 

(226) 

Однородный  сфероид,  касающийся  экватора  Солнца  и  даю- 
щий изопотенциальную  поверхность,  совпадающую  на  эква- 
торе и  полюсе  с  поверхностью  Солнца,  получится  при  кН  =  1. 
Сжатие  этого  сфероида  больше,   чем  у  Солнца: 

_5 

с—  3  ^0- 

Таким  образом,  солнечная  поверхность  будет  внешней  по 
отношению  к  этому  сфероиду,  и  переход  к  случаю  кН  =  1  вполне 
законен,  так  как  выражение  (10)  относилось  к  потенциалу  на 
внешнюю  точку. 

Можно  также  заметить,  что  различие  между  поверхностью 
сфероида  и  его  потенциальной  поверхностью  невелико  и  даже 
на  полюсе  достигает  всего 

.       |аоСо<0",07. 

§  5.  Вращающийся  однородный  сфероид, 

дающий  эквипотенциальную  поверхность  силы  тяжести, 

практически  совпадающую  с  поверхностью  Солнца 

Для  однородного  несжимаемого  сфероида,  вращающегося 
как  твердое  тело,  как  известно,  внешняя  уровенная  поверх- 
ность силы  тяжести  (равнодействующая  силы  тяготения  и 
центробежной  силы)  будет  совпадать  с  его  внешней  поверх- 
ностью, и  эллипсоид  будет  находиться  в  равновесии,  если  между 
его  плотностью  р,  угловой  скоростью  (1>  и  вторьш  эксцентри- 
ситетом 

у'аЗ  —  62 

^=— г- 

2     Вопросы   космогонии,  том   IV 
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существует  соотношение  [61: 

(о2  3  +  12  .  3  /00\ 

2^  =  ^^-агсЬд1-^  =  /..  (23) 


Разлагая   это   выражение  в   ряд,   получим 

00 

4п 


2я^р 


2(-<)-^'4(.  +  1Я-1'"  (23а) 


ИЛИ 

.4  =  Г5'^-Г5'*  +  ^''+--  (236) 

Поскольку  мы  имеем  знакопеременный  сходящийся  ряд, 
то  со  вполне  достаточной  для  нас  точностью  можем  ограни- 
читься в  нем  только  первым  членом  (так  как  1^  ~  10~*).  Тогда 

"^  '^^\  (24) 


27г^р        15 

а  так  как  масса  сфероида 

М  =  4  тгра^б  =  у  тграз  (1  +  Г^)"'^' , 
то 


3   ^Л/         15  21 

или  приближенно 

1-2  _ 

2   ёМ  ' 


'-'^'       '^-'1^  +  ...  (25) 


,•2 5   (О^ДЗ 


Это  уравнение    было   выведено  еще  Ньютоном.'   Вводя  в  (25) 
вместо  второго  эксцентриситета  сжатие  сфероида  с,'  причем 

^■'  =  2с^^^^^^^  =  2с  +  Зс2  +  4сЗ+..., 
получим 

(о2аЗ         4        ,    22    2    ,  ,ог:    ч 

^  =  Т^  +  35^+---  (25а) 

или  приближенно 

_  5^0)203  (26а) 

^  ~"  4  ^Л/  • 
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Если  бы  такой  сфероид  вращался  с  угловой  скоростью, 
равной  угловой  скорости  вращения  поверхностных  слоев 
Солнца 

«0  =  2,868. 10-«, 

то  его  сжатие  было  бы 

с  =  0,2607.  10-*, 

т.  е.  его  сжатие  было  бы  несколько  меньше,  чем  принятое 
нами  наибольшее  сжатие  Солнца  с^. 

Если  мы  возьмем  однородный  сфероид  с  сжатием  Солнца 

Со  =  1.0. 10-*, 

то  по  формуле  (26а)  он  должен  был  бы  иметь  угловую  скорость 
(0  =  1,958(0^  =  5,616.10-8.    . 

Момент  количества  движения  Ь  такого  эллипсоида  был  бы 
значительно  больше  момента  количества  движения  Солнца, 
согласно  модели  Шварцшильда,  вращающегося  как  твердое 
тело  с  угловой  скоростью  (о^:  во-первых,  потому,  что  момент 
инерции  однородного  сфероида 

/  =  |з/а2  (27) 

больше  момента  инерции  модели  Шварцшильда  в  \1д,  или 
в  7,788  раза;  а  во-вторых,  потому,  что  угловая  скорость  его 
больше  почти  в  2  раза.  Момент  количества  движения  получается 
преувеличенным  в  15,23  раза.  Как  мы  видели  в  предыдущем 
параграфе,  такой  сфероид  с  сжатием  Со=1,0-  10-*  будет  иметь 
уровенную  поверхность  силы  тяготения  несколько  менее  сжа- 
тую, чем  его  собственная  поверхность. 

Мы  желаем,  чтобы  потенциальная  поверхность  силы  тяго- 
тения имела  то  же  сжатие,  что  и  внешняя  поверхность  Солнца. 
Соответствующий  однородный  сфероид,  как  мы  видели  в  пре- 
дыдущем параграфе,  должен  иметь  сжатие,  равное 

с  =  |со  =  1,667.10-*. 

Если  мы  в  формулу  (26а)  подставим  это  значение  с,  то  получим 
(5ще  большую  угловую  скорость,  равную 

О)  =  2,529(00=  7,252.  10-е. 
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Соответствующее  преувеличение  момента  количества  движения 
будет  в  19,70  раза: 

^^  =  19,70^о. 

§  6.  Модель  с  однородным  ядром  и  невесомой 

оболочкой  » 

По  известной  теореме  Маклорена,  однородные  софокусные 
эллипсоиды  с  совпадающими  центрами  и  совпадающим  направле- 
нием осей  имеют  одинаковые  внешние  уровенные  поверхности 
потенциала  тяготения.  Если  их  массы  равны,  то  общей  будет 
не  только  форма  уроненной  поверхности,  но  и  значение  потен- 
циала на  ней.  Мы  рассмотрим  в  этом  параграфе  семейство 
моделей  Солнца,  состоящих  из  вращающегося  ядра  (разных 
размеров),  содержащего  всю  массу  Солнца,  и  из  невесомой 
оболочки.  Ядро  будем  представлять  себе  однородным  и  софо- 
кусным  с  эллипсоидом,  касающимся  экватора  Солнца  и  имею- 
щим сжатие 

5 

Такие  эллипсоиды  будут  давать  внешнюю  потенциальную 
поверхность  силы  тяготения,  практически  совпадающую  с  по- 
верхностью Солнца  и  имеющую  то  же  сжатие.  Сама  поверх- 
ность наибольшего  эллипсоида  практически  также  будет  сов- 
падать с  поверхностью  Солнца,  так  как  его  полярная  полуось  Ь^ 
будет  отличаться  от  полярной  полуоси  Солнца  \  только  на  О", 07, 
Для  софокусных  сфероидов 

а2  — 62  =  а2  — 62 
и 

е2^е2-^  =  ^-  — --^  Г28) 

где  эксцентриситет  наибольшего  сфероида 

2         1,2 

„о О  1 

е,  =   г 


И,  как  и  прежде, 


а^к. 


формула  (226)  §  4  является,  таким  образом,  следствием  теоремы 
Маклорена,  если  мы  учтем,  что 


3    о- 


0  моменте  количества  движения 

Солнца 
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Значения  с,  е^  и  1^  для  софокусных  сфероидов  при  различных 

значениях  А:  приведены  в  табл.  3. 

Таб 

лица    3 

С:исатие  сфероидов,    имеющих    эквипотенциальную    поверхность,    совпа- 

дающую с  поверхностью  Солнца 

«0 

с- 10* 

е'  •  10* 

г' .  10< 

Ш/Шд 

ыи 

1,0 

1,667 

3,334 

3,335 

2,5 

1,000 

0,9 

2,058 

4,116 

4,118 

3,5 

1,054 

0,8 

2,605 

5,210 

5,213 

4,4 

1,117 

0,7 

3,402 

6,804 

6,809 

6,2 

1,192 

0,6 

4,631 

9,262 

9,271 

9,1 

1,274 

0,5 

6,668 

13,336 

13,354 

14,3 

1,345 

0,4 

10,419 

20,838 

20,882 

25,0 

1,358 

0,3 

18,522 

37,044 

37,182 

51,3 

1,184 

0,2 

41,675 

83,350 

84,051 

149,4 

0,760 

0,12 

115,76 

231,52 

237,01 

506,8 

0,454 

Как  видно,  допустимое  сжатие  этих  сфероидов  с^ень  мало, 
и  даже  для  случая  всей  массы,  заключенной  в  объеме  конвек- 
тивного   ядра    Солнца,    сжатие 

с<0,0И6. 
Момент  количества  движения  таких  эллипсоидов  равен 

Ь  =  /(О , 

где  /  —  момент   инерции  —  дается  формулой  (27),   а   угловая 
скорость  О)  может  быть  определена  из  формулы  Ньютона  (26а): 


4  ^М 

5  "^ 


^.=... 


Более  точно  м  может  быть  найдено  из  формулы  (25а) 
Из  (29)  и  (22)  получим 


(29) 


(29а) 


'-~У    3   оЗ    А2 
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н,  так  как 


О' 


„=(|«|У'^-о;\  (30) 

то  момент  количества   движения  для    рассматриваемых   моде- 
лей будет 

Согласно  (27),  мы  можем  считать 

*  5         0 


^=^V^м\;ч-'^  (31) 

При  Л=1  мы  получим 

Величины  ЫЬ^  приведены  в  последнем  столбце  табл.  3, 
а  значенЕЙ  м/мо,  вычисленные  по  формуле  (30),  —  в  пред- 
последнем. Как  видно  из  этой  таблицы,  увеличение  концен- 
трации к  центру  сопровождается  очень  большим  увеличением 
угловой  скорости,  но  сравнительно  небольшим  увеличением 
момента  количества  движения.  Наибольшая  скорость  вращения 
дает  конфигурацию,  еще  далекую  от  формы,  переходной  от 
эллипсоидов  Маклорена  к  эллипсоидам  Якоби. 

§  7.  Двухслойная  модель  Солнца 

В  предыдущем  параграфе  мы  рассматривали  Солнце  как 
состоящее  из  однородного  сжатого  ядра,  быстро  вращающегося, 
и  невесомой  оболочки.  Считая  ядро  однородным,  а  не  уплот- 
ненным к  центру  и  сосредоточивая  в  нем  всю  массу  Солнца, 
мы  преувеличивали  момент  инерции,  а  следовательно,  и  мо- 
мент количества  движения.  Но,  пренебрегая  массой  оболочки 
и  ее  вращением,  мы  несколько  уменьшали  момент  количества 
движения.  Если  быстро  вращающиеся  части  Солнца  имеют 
радиус,  больший  0,5  радиуса  Солнца,  то  мы,  повидимому, 
получаем  преувеличенное  значение  для  момента  количества 
движения,  поскольку  масса  частей  Солнца  вне  сферы  с  /-=0,5  г- 
составляет   менее   4%    всей   массы   Солнца.    При   значительно 
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меньших  размерах  вращающегося  ядра  наша  приближенная 
оценка  «сверху»  момента  количества  движения  не  будет  уже 
годиться. 

Рассмотрим  теперь  двухслойную  модель  с  граничной  поверх- 
ностью в  форме  того  же  софокусного  эллипсоида.  Внепший 
слой  будем  считать  вращающимся  с  угловой  скоростью  поверх- 
ностных слоев  (Од,  а  внутреннее  ядро  —  с  угловой  скоростью, 
определяемой  по  формуле  Ньютона  (26а),  хотя  эта  формула 
здесь  и  не  будет  строго  применима,  поскольку  поверхность 
ядра  не  является  свободной  поверхностью,  а  на  нее  влияет 
потенциал  внешних  слоев. 

Пусть  масса  эллипсоида  с  большой  полуосью 

будет  та  же,  как  и  для  сферы  того  же  объема  в  модели  Шварц- 
шильда.  Тогда 

г;  =  а^(1-с)  =  аг^(1-|-|-);         '  (32) 

г,  =  а/(1-|^)''«аоЛ-(1-|-^);  (32а) 

о 

Величина  М^  может  быть  взята  из  табл.  1. 
Средняя  плотность  центрального  ядра  будет 

Рв=^Л/(г.).  (33) 

Средняя  плотность  оболочки  будет  равна 

_  М  (гр)  —  М  (Г1)  .  .о/ч 

момент  инерции  оболочки 

Iи=^[М{^о)-^Ц^,)]а^{^-к^^).  (35) 


24  Н.  Н.  Парийский 


Наконец,  момент  инерции  ядра  будет  равен 

1,  =  ^М{г,)а1Е\  (36) 

момент  количества  движения  двухслойной  модели  выразится 
суммой 

^  =  /д  0)0  -1-  /в  <о ,  (37) 

где  «о  —  угловая  скорость  поверхностных  слоев,  а  со  —  угло- 
вая скорость  ядра,  определяемая  по  формуле  Ньютона  через 
сжатие  граничного  сфероида. 

Если  этот  граничный  сфероид  является  софокусным  с  наи- 
большим однородным  сфероидом,  определенным  в  предыдущем 
параграфе,  то  внешнее  поле  от  этих  двух  сфероидов  будет 
таким  же,  как  и  от  одного  наибольшего,  и  даст  такое  же  сжа- 
тие внешней      уроненной  поверхности. 

Действительно,  мы  можем  применить  теорему  Маклорена 
о  внешнем  поле  однородных  сфероидов,  рассматривая  нашу 
модель  как  случай  однородных  сфероидов,  наложенных  друг 
на  друга,  с  плотностями  рн  и  рв  —  рн  и  массами,  в  сумме 
равными  массе  всего  Солнца. 

Угловая  скорость  со,  определяемая  по  формуле  (30),  в  этой 
модели  будет  меньше,  чем  в  предыдуш,ей,  так  как  под  М  мы 
будем  теперь  подразумевать  только  массу,  заключенную  внутри 
эллипсоида : 


=--)/ 


4      Л/ в 


"  =  \/|  ^  \1с,  к-"'  \1М,  ■=  со,  |/^  =  со,  фМ, ,         (38) 

м 
где  (01  дается  формулой  (30),  а  М^=^  приведено    в  табл.  1. 

В  табл.  4  приведены  результаты  вычислений  момента  коли- 
чества движения  Солнца  для  различных  положений  границы 
раздела,  для  различных  значений  к.  В  первом  столбце  даны 
значения  к=а1а^,  во  втором  —  радиусы  сфер,  вычисленные  по 
формуле  (32)  в  единицах  радиуса  Солнца;  в  третьем  —  массы 
внутренних  эллипсоидов  в  долях  массы  Солнца.  Обозначение 
последующих  столбцов  ясно  из  текста. 

Как  видно  из  табл.  4,  момент  количества  движения  Солнца, 
согласно  этой  модели,  также  не  намного  может  превосхо- 
дить момент  для  Солнца,  вращающегося  как  твердое  тело. 
Любопытно,  что  наибольший  момент  количества  движения  из 
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Таблица   4 
Момент  количества  движения  Солнца  для  двухслойной  модели 


к 

М(г,) 

М(Го) 

Рн 

/в 

/и  •  10~53 

/3-10-53 

ш.]0' 

(и 

«"о 

Ь 

1,00 

0,99994 

1,0000 

0,0000 

1,4112 

0,000 

38,333 

7,250 

2,528 

19,69 

0,90 

0,89994 

1,0000 

0,0000 

1,9358 

0,000 

31,050 

9,437 

3,291 

20,76 

0,80 

0,79993 

0,9997 

0,0009 

2,7555 

0,004 

24,525 

12,660 

4,415 

21,99 

0,70 

0,69992 

0,998 

0,0043 

4,1067 

0,039 

18,745 

17,664 

6,160 

23,46 

0,60 

0,59991 

0,990 

0,0180 

0,4697 

0,245 

13,662 

25,861 

9,018 

25,08 

0,50 

0,49989 

0,962 

0,0613 

10,865 

1,092 

9,219 

40,218 

14,025 

26,49 

0,40 

0,39986 

0,887 

0,1703 

19,576 

3,638 

5,439 

67,467 

23,527 

26,73 

0,30 

0,29981 

0,706 

0,4263 

36,961 

10,273 

2,434 

123,102 

42,929 

23,31 

0,20 

0,19957 

0,380 

0,8818 

67,440 

22,808 

0,583 

249,872 

87,136 

14,95 

0,12 

0,11954 

0,112 

1,2550 

92,410 

33,541 

0,617 

486,382 

169,61 

8,94 

этой  серии  моделей  получается  для  модели,  в  которой  цен- 
тральное ядро  занимает  0,4  радиуса  Солнца.  В  этом  случае 
ядро  содержит  89 °о  Бсей  массы,  97,2%  всего  момента  коли- 
чества движения  и  вращается  с  угловой  скоростью,  в  9,17  раза 
превышающей  угловую  скорость  оболочки.  Момент  количества 
движения  в  этой  модели  «наибольшего  момента»  в  26,7  раза 
больше,  чем  для  первой  из  рассмотренных  нами  моделей,  при- 
чем, как  уже  указывалось,  основное  увеличение  происходит 
вследствие  принятой  однородности  оболочки  и  ядра  и  вслед- 
ствие увеличения  угловой  скорости. 

§  8.  Случай  дифференциального  вращения 
многослойной  модели  Солнца 

В  предыдущем  параграфе  мы  рассмотрели  двухслойную  мо- 
дель Солнца.  Она  может  дать  только  очень  грубую  оценку,  — 
невидимому ,  с  некоторым  преувеличением  допустимого  момента 
количества  движения  Солнца. 

Можно  рассмотреть  и  такую  модель  Солнца  с  непрерывно 
меняющейся  плотностью,  чтобы  внешняя  уровенная  поверх- 
ность совпадала  с  поверхностью  Солнца  на  полюсе  и  экваторе. 
Для  этого  представим  себе  Со;у1це  состоящим  из  бесконечного 
числа  бесконечно  тонких  слоев,  заключенных  между  софокус- 
ными  эллипсоидальными  поверхностями,  причем  каждый  слой 
будем  считать  вращающимся  со  своей  угловой  скоростью,  опре- 
деляемой по  формуле  Ньютона  (26а).  Этим  мы  вводим  несколько 
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произвольный  закон  распределения  скоростей,  поскольку  фор- 
мула Ньютона  не  строго  применима  к  этой  модели.  Однако 
можно  ожидать,  что  с  точностью  до  порядка  мы  все  же  полу- 
чим оценку  допустимого  момента  количества  движения. 

Такую  модель  можно  также  представить  себе  состоящей  из 
наложения      бесконечного      числа      однородных     эллипсоидов 

с  плотностью   "о"  ^  с?Г1 ,   где   г^  определяется    уравнением    (32) 

или  (32а). 

Представим  себе  сперва  стандартную  сферическую  модель 
внутреннего  строения  Солнца,  например  модель  Шварцшйльда. 
Распределение  всех  интересующих  нас  функций  в  этой  модели 
приведено  в  табл.  1.  Деформируем  теперь  сферические  слои 
этой  модели  так,  чтобы  слой,  заключенный  между  сферами 
с  радиусами  г  и  г-\-Лг,  превратился  в  сфероидальный  слой 
между  двумя  софокусными  эллипсоидами  с  тем  же  объемом 
и  с  теми  же  массой  и  плотностью.  Внешняя  потенциальная 
поверхность  такой  деформированной  сферы  будет  опреде- 
ляться деформацией  наружной  поверхности  сферы,  так  как 
эллипсоиды  софокусные. 

Мы  будем  считать  потенциальную  поверхность,  касающуюся 

экватора  деформированной  сферы,   имеющей  сжатие  с^.  Тогда 

наружная  эллипсоидальная  поверхность  будет  иметь,  согласно 

5 
(22),    сжатие    с^  =  ^Со,    а    большие   полуоси  соответствующих 

эллипсоидов  будут  определяться  из  уравнений  (32)  или  (32а), 
или  а=ка^,  где  к  определяется  из  уравнения: 

^(1-||^)=^^-'  (39) 

"о 

Наибольшее  значение  г  лишь  немного  меньше  а^: 

|=(1-А.,).  (39') 


Сжатия  софокусрых  эллипсоидов  будем  считать  приблизи- 
тельно поддерживающимися  их  вращением.  Рассмотрим  один 
из  внутренних  эллипсоидов.  Внутри  него  находится  масса 
М{г).  Условие,  что  мы  считаем  эту  массу  распределенной 
по  эллипсоидальным  софокусн&м  слоям,  не  влияет  на  внешнее 
поле  этого  эллипсоида.  Точно  так  же  и  условие,  что  отдельные 
слои  вращаются,  тоже  не  влияет  на  потенциал  на  поверхности 
этого  слоя.  На  потенциал  в  слое  имеет  влияние  только  угловая 
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скорость  вращения  самого  слоя.  Мы  видели,  что  уровенная 
поверхность  силы  тяготения  такого  эллипсоида  будет  иметь 
сжатие,  меньшее,  чем  сам  эллипсоид  (именно,  ^5  от  его  сжатия). 
Мы  считаем,  что  уровенная  поверхность  совпадает  с  поверх- 
ностью эллипсоида.  При  этом  мы  делаем  ошибку  порядка 
величины  сжатия  эллипсоида.  Можно  заметить,  что  учет  потен- 
циала центробежной  силы  будет  делать  сжатие  внешней  уро- 
ненной поверхности  большим,  а  учет  потенциала  наружных 
слоев  будет  его  уменьшать.  Таким  образом,  можно  думать, 
что  рассматриваемая  модель  не  спльно  отличается  от  динами- 
чески возможной.  Нахождение  динамической  модели  равно- 
весной конфигурации  с  дифференциальным  вращением  слоев 
представляет  задачу,  насколько  нам  известно,  еще  не  решен- 
ную, если  не  считать  работ  Дирихле  и  Римана  для  гомологи- 
ческих движений  в  однородных  эллипсоидах.  Поскольку  нас 
интересует  только  приблизительный  подсчет  возможного  мо- 
мента количества  движения  Ск)лнца,  мы  ограничимся  нашей 
моделью,  которая,  как  нам  представляется,  не  должна  давать 
момент  количества  движения,  отличающийся  во  много  раз  от 
динамически  равновесной  модели. 

Уравнение  изопотенциальной  поверхности  потенциала  тяго- 
тения для  модели  из  вложенных  друг  в  друга  однородных  софо- 
кусных  сфероидов  может  быть  написано  в  таком  виде: 

*"© 
У{х,  у,  г)^о{х,  у,1)§\  аМ,  (40) 

о 

где  о  {х,  у,  2)  —  функция,  определяющая  форму  у ро венной 
поверхности,  одинаковую  для  каждого  из  вложенных  друг 
в  друга  однородных  сфероидов  (в  дальнейшем  мы  будем  назы- 
вать эти  сфероиды  элементарными);  (1М  —  масса  элементар- 
ного сфероида. 

Размерность  функции  о  {х,  у,  г)  будет  см-^ 
Масса    элементарного    сфероида    равна    произведению    его 
плотности  р  на  объем  V: 

где,  очевидно,  по  условию  равенства  объемов  сфероида  и  соот- 
ветствующей сферы. 

Плотность  элементарного  сфероида,  очевидно,  должна  быть 
взята  равной  средней  избыточной  плотности  между  соседними 
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сфероидальными  поверхностями  или,  по  условию  задачи, 
между  соседними  сферическими  слоями,  т.  е. 

Действительно,  — ^(1г  есть  разность  плотностей  между  сосед- 
ними слоями  при  перемещении  снаружи  внутрь  Солнца,  и  сред- 
няя избыточная  плотность  должна  равняться  половине  этой 
величины.   Таким  образом, 

ам=-^^г^'/^ап  (42) 

У{х,  у,  г)  =  -о(х,  у,  г)  \   ^-г'^^с1г.  (43) 

о 

Нас  интересует  момент  количества  движения  для  нашей 
модели.  Этот  момент  количества  движения  может  быть  пред- 
ставлен как  сумма  моментов  количества  движения  отдельных 
элементарных  сфероидальных  слоев: 

г 

о 

где   1{г)йг  —  момент   инерции    элементарного  слоя,   а  (о  (г)  — 
его  угловая  скорость. 
Легко  видеть,  что 

/  (г)  аг  =  -|  р  (г)  \~т:{г-\-  йгУ  а^  (г  +  а,')  -  ^  ^г^а^  (г)]  , 
/(г)с/г=^^гр(г)г^[(1  +  ^^/а2(г  +  ^г)-аМг)]. 


или 


Но 

а  (г)  =  «/(/•); 


аЦг-^аг)  =  а1кЦг-^с1г)^а1\_к{г)-^''-^Лг'^== 

/(^)^^=^.р(.)а^А-(.)^^[(1  +  ^^)(1+247^^0-^]' 

/  (г)  Лг=  ^.  .0^  (г)  гЧ^-  (г)  [1  +  ^  ^]  аг.  (44) 
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с1к{г)1йг  может  быть   найдено  из  уравнения  (39),   определяю- 
щего к{т).  Из  этого  уравнения  получаем 

Ограничиваясь  двумя  членами  и  дифференцируя  (45), 
получаем 

^  _.  _1_  Л    I    _5  со^ 12  ^^\ 

В  малых  членах  правой  части  заменяем 

к  =  ^- 
ч  ' 
тогда 

^  =  ^['-|^оШ         '  .(46) 

Из  (45)  с  достаточным  приближением  имеем 

Подставляя  это  значение  к{г)  в  (44),  получим 

/  (г)  с//-  ==  1  гр  (г)  г*  [1+  I  Со (^)']  ^г. 

Легко  видеть,  что  тг  -о  {г)  г^Лг  есть  момент  инерции  сфе- 
рического слоя  с  плотностью  р(г),  радиусом  г  и  толщиной  йг, 
а      1-[--^Со(— ]     — поправочный    множитель   для   получения 

момента  инерции  сфероидального  слоя  между  двумя  софокус- 
ными  сфероидальными  поверхностями. 

Угловую  скорость  вращения  элементарного  слоя  мы  опять 
примем  условно  определяемой  уравнением  (30),  где  положим 

М  =  М{г) 
—  массе,  заключенной  внутри  радиуса  г;  тогда  из  (30)  и  (47) 


или 


('•)  =  (I  «)''■  %М  М'Ч .г--'.  [1-1  с„(^У ]  ,  (48) 
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ИЛИ 


(Н^»Г*^Ш'"1^-^-^^]-  (^««) 


('•)=  т-1 


Момент   количества  движения  элементарного  сфероидального 
слоя  будет 

/  (г)  (О  (г)  (1г  = 

=  ^  7г^'/.а^с;/^  [Ж  (,)]-/.р  (,)  г^и  [1  _  1|  ,^ (^у]  ,/,,        (49) 

а  полный  момент  количества  движения  этой  модели  будет 

'о 


Вводя 
получим 


•=''1''0 


Лиг 


13  Со"!  л 


(50) 


О 

Полагая 

с©  ==1,0.  10-* 

и  пользуясь  данными  табл.  1  для  М {г^)  и  р(^1),  мы  получаем 
^  =  6,108. 10*8  смЧек-!. 

Таким  образом,  для  этой  модели,  более  близкой  к  равно- 
весным динамическим  моделям,  момент  количества  движения 
Солнца  получается  всего  в  4,33  раза  больше,  чем  для  Солнца, 
согласно  модели  Шварцшильда,  вращающегося  как  твердое 
тело  с  угловой  скоростью  ыо- 

Распределение  момента  количества  движения  Солнца  в  этой 
модели  по  отдельным  слоям,  определяемое  формулой  (49),  при- 
ведено в  табл.  5,  во  втором  столбце.  В  третьем  столбце  приве- 
дены угловые  скорости  {о/сод  по  отдельным  слоям. 

В  последнем  столбце  приведены  градиенты  скорости  Л'в/б?г 
в  км/сек  на  1  км,  где  гв  —  круговая  скорость  в  экваториальной 
плоскости,  или,  что  то  же,  в  см/сек  на  1  см.  Очевидно,  что 


а»0 аыг ^^     010 

(1г  (1г  '     '       с1г 


(51) 
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Таблица  5 

Распределение      момента     количества    двимсения,      угловых 

скоростей    и    градиентов    скоростей    по    отдельным    слоям 

в  модели    с  дифференциальным  вращением 


г 

I  (г)  «о  (г)  йг  •  10-*в 

»(г) 

йг 

0,00 

0 

со 

00 

0,10 

368 

208,4 

-6,1 

0,15 

869 

129,2 

—10,2 

0,20 

1187 

86,8 

—12,8 

0,25 

1152 

60,4 

—13,0 

0,30 

906 

43,0 

—11,8 

0,35 

627 

31,4 

—10,0 

0,40 

401 

•    23,5 

-8,3 

0,45 

243 

18,0 

—6,8 

0,50 

142 

14,0 

—5,6 

0,55 

80 

11,2 

-4,5 

0,60 

43 

9,0 

-3,8 

0,65 

22 

7,4 

-3,2 

0,70 

10 

6,2 

-2,5 

0,75 

5 

5,2 

-2,3 

0,80 

1,7 

4,4 

-1,9 

0,85 

0,5 

3,8 

-1,7 

0,90 

ОД 

2,9 

-1,5 

0,95 

0,02 

2,9 

-1,2 

1,00 

0,00 

2,5 

-1,0 

Как  видно,  эти  градиенты  очень  малы.  Вопроса  о  влиянии 
вязкости  и,  в  частности,  лучистой  вязкости  на  замедление  диф- 
ференциального вращения  мы  в  данной  работе  не  касаемся. 

Легко  вычислить  для  последней  модели  механическое  сжатие 
Солнца, 

Момент  инерции  относительно  экваториальной  оси  для  такой 
модели  будет,  очевидно, 


Но 
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а  е^  определяется  из  уравнения  (28): 

о2. 


ъ  к^ 


кроме  того, 


а''  =  к^{г)а1^ 


причем   А:^   определяется  из  уравнения  (39).    Таким   образом, 

О 

Подставляя  значение  йМ  из  (42),  получим 

Д  =  +  А  ,„.  |\-.  (,)  [1  _  1  _|2_]  ^1  ^1Р  I  ^,.  (52) 

О 

Этот  момент  инерции  будет  мало  отличаться  от  момента  инер- 
ции модели  Шварцшильда. 

Момент  инерции  относительно  полярной  оси  будет 


I  а^ам 


или 


••© 


/р=Г5^«?1   ^"^')'' 


аг 


(53) 


Определяем  динамическое  сжатие  для  нашей  модели: 


^Р-^. 


6    ]         \с1г 


дин 


I     кЦг)г 


аг 


аг 


(54) 


Вычисляя  (1о/с1г  численным  методом  по  данным  табл.  1  и 
графическим  методом  для  концов  интервала,  можно  опреде- 
лить 


/„  =  2,454.10 
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Таким  образом,  механическое  сжатие  для  этой  модели 
всего  в  6,5  раза  больше,  нежели  сжатие  внешней  эквипотен- 
циальной поверхности  силы  тяготения.  Произведенные  расчеты 
показывают,  что  момент  количества  движения  Солнца,  даже  при 
допущении  существования  быстро  вращающихся  внутренних 
слоев,  может  составлять  только  около  2%  момента  количества 
движения  всех  планет  солнечной  системы. 

Усложнение  рассмотренной  нами  модели  вряд  ли  сможет 
изменить  эту  цифру  более  чем  в  2  раза.  Таким  образом,  воз- 
можность падения  на  Солнце  и  стягивания  в  Солнце  сколько- 
нибудь  значительной  доли  материала,  послужившего  для  обра- 
зования планет,  если  не  учитывать  потерь  Солнцем  своего 
момента  количества  движения  каким-либо  путем  и  если  счи- 
тать, что  этот  материал  вращался  вокруг  Солнца  подобно 
планетам,  должна  быть  исключена. 
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